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ABSTRACT 


This  report  presents  an  algorithm  and  HP-41CV  programs  for 
position  determination  with  Loran-C  chains.  Additional 
computational  routines  include  the  ability  to  calibrate  Loran 
station  triplet  data  to  a  known  benchmark  and  ITD's  (Indicated 
Time  Delay's),  predict  ITD's  at  given  positions,  compute  the 
geodesic  (similiar  to  great  circle)  bearing  and  distance  from  a 
fix  to  the  destination  and  to  compute  the  geodesic  bearing  and 
distance  from  any  one  location  to  another.  Utility  routines  allow 
the  user  to  transfer  station  pair  data  between  the  HP-41CV  and 
magnetic  cards,  magnetic  tape  and  an  extended  function/  memory 
module. 


a.   introduction 

The  Loran  system  is  a  radio  aid  to  navigation  which  utilizes 
the  principle  of  hyperbolic  fixing.  The  locus  of  points  for  which 
the  difference  in  arrival  time  of  synchronized  signals  from  a  pair 
of  transmitters  is  constant  determines  a  hyperbolic  line  of 
positions.  The  intersection  of  two  hyperbolic  lines  of  position 
from  two  pairs  of  stations  determines  position  or  a  hyperbolic 
fix.  That  two  pairs  of  stations  are  required  for  a  fix  does  not 
necessarily  mean  that  there  are  four  separate  stations,  for  one 
station  of  one  pair  may  be  colocated  with  one  station  of  the  other 
pair  forming  a  Loran  triplet  (Figure  1).  Triplets  may  be  joined 
"end-to-end"  by  station  colocation  to  form  a  Loran  chain  (Figure 
2).  Loran  chains  are  common  on  both  the  East  and  West  coasts  of 
the  North  American  continent. 

The  present  day  Loran-C  operates  at  100-kHz  and  is  in  use  in 
the  Atlantic,  Pacific  and  Mediterranean  areas.  The  computational 
algorithm  and  programs  described  herein  can  be  used  for  position 
determination  with  Loran-C  triplets.  Further  information  on  the 
history,  development  and  operation  of  the  Loran  systems  may  be 
found  in  References  1  and  2. 


MrM2 


(a)   Colocated  Master  Stations 


(b)   Colocated  Slave  Stations 
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(c)   Colocated  Master  and  Slave 


Figure  1.   Loran  Triplets 
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Figure  2.   Loran  Chain  of  Five  Loran  Triplets 


b.  us£L  Routines 

1.  Computational  Routines: 

FI  -  The  F_Ixing  routine  is  the  main  program  for  calculating 
a  Loran-C  fix  from  indicated  time  delays. 

AS  -  The  Alternate  Solution  routine  will  allow  the  second 

Loran  fix  solution  to  be  computed.  This  routine  toggles 
Flag  3  so  that  on  subsequent  fixes  the  FI  routine  will 
calculate  the  alternate  solution. 

DJi  -  The  Destination  routine  stores  the  latitude  and 
longitude  of  a  fixed  destination. 

HD_  -  Computes  the  Heading  and  distance  from  the  current  fix 
to  the  destination  stored  by  the  DN  routine. 

2.  Manual  Mode  Routines: 

Ml  -  Manual  Input  allows  station  data  to  be  input  and  stored 

via  the  calculator  keyboard. 
ED  -  £cho  D_ata  is  a  utility  routine  for  validating  station 

triplet  information  stored  in  the  calculator. 

3.  Card  Reader  Routines: 

££  -  £ard  Store  records  station  data  onto  magnetic  cards. 
QR  -   £ard  Bead  inputs  station  data  from  magnetic  cards. 
£_E  -  £ard  .Echo  is  a  utility  routine  for  validating  station 
information  stored  on  data  cards. 

4.  Extended  Memory  Routines: 

X£  -  XMEM  Store  records  station  data  onto  the  extended 

memory  module. 
XR  -   XMEM  Becall  inputs  station  data  from  the  extended 

memory  module. 


XD  -  XMEM  Delete  erases  station  data  from  the  extended 
memory  module. 

5.  Tape  Cassette  Routines: 

T£  -  Tape  S_tore  records  station  data  onto  the  tape  cassette. 
IB  -  T_ape  Becall  inputs  station  data  from  the  tape  cassette. 
ID.  -   Tape  Delete  erases  station  data  from  the  tape  cassette. 

6.  Utility  Routines: 

DH  -  This  routine  is  similiar  to  the  HD  option  except  that 

it  computes  the  heading  and  distance  from  any  origin  to 
any  destination. 

£E  -  PR  is  used  to  Predict  the  station  ITD's  that  will  be 
received  at  a  given  latitude  and  longitude. 

£A.  -  CA  is  the  Calibration  option.  Given  the  latitude  and 
longitude  of  a  known  position  and  the  indicated  time 
delays  from  a  Loran-C  triplet,  the  stored  station  data 
are  modified  so  that  the  FI  routine  (or  AS)  will 
compute  the  known  position  from  the  same  time  delays. 

SW  -  .Sjtfitch  data  swaps  the  data  of  the  two  Loran  stations 
stored  in  memory. 

Note:  There  are  no  specific  routines  that  relate  to  HP-41CV 
printer  operations.  However,  all  input  and  output  will  be 
recorded  on  a  printer  if  one  is  attached. 


c.  Recommended  Key  Assignments 

It  is  recommended  that  the  following  HP-41CV  Ais_e_r_  key 
assignments  be  made  and  recorded  onto  the  program  cards  if  the 
program  is  first  prepared  manually: 


CA   ED 


cs 

CR 

CE 

TD 

MI 

sw 

TS 

TR 

DN 

XD 

PR 

DH 

xs 

XR 

HD 

A 

S 

FI 

The  £1  function  is  placed  on  the  ENTER  key  as  a  reminder  that 
ENTER  should  not  be  used  for  data  entry.  To  make  the  user  key 
assignments  on  the  HP-41CV,  refer  to  the  ASN  function  in  the 
HP-41CV  Owner's  Handbook  and  Programming  Guide. 


D.  User  Instructions  with  Examples 
1.  Manual  Input  Routine  -  MI 

This  routine  can  be  used  to  enter  station  data  (Appendix 
C)  to  prepare  data  to  be  transferred  to  cards,  extended  memory, 
or  tape  cassette.  It  can  also  be  used  to  enter  station  data 
manually  if  alternate  data  storage  media  are  not  available. 

Example  using  9940W  and  9940Y.  (Note:  The  notation  <9>, 
for  example,  means  that  you  must  press  the  gold  shift  key  and 
tJie_n  the  9  key.) 

£ej£  Key  in        Rlsss. 

1.  XEQnMI" 

2.  ID?  <9><9><4><0>W  R/S 

3.  CODE  DELAY?         11000  R/S 

4.  MS  LAT?  39.330662  R/S 

5.  MS  LON?  118.495637  R/S 

6.  SS  LAT?  47.034799  R/S 

7.  SS  LON?  119.443953  R/S 

If  desired,  these  station  parameters  can  now  be 
transferred  to  card,  extended  memory  or  tape  using  the  CS,  XS, 
or  TS  routines,  respectively.  Otherwise,  repeat  the  steps  above 
with  the  9940Y  station  data.  These  two  stations  will  be  used  in 
the  remaining  examples. 

£££         Key  in        £te_s_s_ 

1.  XEQ"MI" 

2.  ID?  <9><9><4><0>Y       R/S 

3.  CODE  DELAY?         40000  R/S 

4.  MS  LAT?  39.330662         R/S 


118.495637 

R/S 

35.191818 

R/S 

114.481743 

R/S 

5.  MS  LON? 

6.  SS  LAT? 

7.  SS  LON? 

8.  NEXT  OPTION? 

At  this  point  the  station  data  for  9940W  and  9940Y  are 
stored  in  the  calculator.  The  9940Y  data  can  be  transferred  to 
card,  extended  memory  or  tape  cassette  using  the  CS,  XS  or  TS 
routines,  respectively. 

[Advanced  User  Note:  The  first  action  taken  in  the  MI 
routine  is  to  XEQ"SWn.  SW  is  the  data  swap  routine  which 
exchanges  to  content  of  R16  -  R25  with  R26  -  R35.  The  incoming 
data  are  then  stored  in  the  R16  -  R25  registers.  The  CS,  XS  and 
TS  routines  transfer  the  content  of  R16  -  R25  to  card,  extended 
memory  or  tape  cassette,  respectively.  If  desired,  the  content 
of  these  registers  can  be  swapped  once  more  by  using  the  SW 
utility  routine.] 


2.  Echo  Data  Routine  -  ED 

This  routine  allows  the  user  to  review  the  station  data 
resident  in  the  calculator. 

Example:  Load  the  station  pairs  9940W  and  9940Y  using 
either  the  MI,  CR,  XR  or  TR  routines. 

See  Key  in         Press 

1.  XEQnED" 

2.  ID:  9940W  R/S 

3.  CD:  11000  R/S 

4.  MLT:  39.330662  R/S 

5.  MLN:  118.495637  R/S 

6.  SLT:  47.034799  R/S 

7.  SLN:  119.443953  R/S 

8.  BL:  2796.903  R/S 

9.  ID:  9940Y  R/S 

10.  CD:  40000  R/S 

11.  MLT:  39.330662  R/S 

12.  MLN:  118.495637  R/S 

13.  SLT:  35.191818  R/S 

14.  SLN:  114.481743  R/S 

15.  BL:  1967.302  R/S 

16.  NEXT  OPTION? 

Notation:  CD  =  coding  delay,  M  =  master,  S  =  slave,  LT  = 
latitude,  LN  =  longitude  and  BL  is  the  station  pair  baseline 
plus  the  secondary  phase  correction  in  microseconds. 


3a.   Card  Store  Routine  -  CS 

With  the  card  reader  attached,  station  data  (in  R16  - 
R25)f  which  has  been  input  using  the  MI,  XR  or  TR  routine,  can 
be  transferred  to  magnetic  card  using  the  XEQ"CS"  command  or  by 
pressing  the  appropriate  user  defined  key. 

Example  using  9940W. 

£££  Key  in  Px£S_s_ 

1.  XEQ"PD" 

2.  WRITE:  9940W    (Pass  a  blank  card  through  the  card 

reader.  Label  the  card  track  n9940Wn) 

3.  NEXT  OPTION? 

To  proceed  with  the  remaining  examples  it  is  recommended 
that  you  also  prepare  a  card  for  the  9940Y  station  pair. 

3b.  Card  Read  Routine  -  CR 

With  the  card  reader  attached,  XEQ"CR"  or  press  the 
appropriate  user  defined  key.  This  routine  can  be  used  to  input 
the  data  for  two  station  pairs,  which  must  form  a  triplet. 

Example  using  9940W  and  9940Y. 

£££  Key  in        p_rje_££ 

1.  XEQnCR" 


2.  1ST  CARD     (Pass  the  data  card  for  9940W  through 

the  card  reader.) 

3.  2ND  CARD     (Pass  the  data  card  for  99440Y  through 

the  card  reader.) 

4.  NEXT  OPTION? 
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3c.   Card  Echo  Routine  -  CE 

With  the  card  reader  attached,  XEQ"CE"  or  press  the 
appropriate  user  defined  key.  This  routine  is  used  to  validate 
the  content  of  data  cards  against  the  table  in  Appendix  C. 

Example  using  9940W. 

£££  Key  in         Press 

1.  xeq-ce" 

2.  STA.  CARD    (Pass  one  side  of  a  data  card  [9940W] 

through  the  card  reader) . 

3.  ID:  9940W  R/S 

4.  CD:  11000  R/S 

5.  MLT:  39.330662  R/S 

6.  MLN:  118.495637  R/S 

7.  SLT:  47.034799  R/S 

8.  SLN:  119.443953  R/S 

9.  BL:  2796.903  R/S 
10.  NEXT  OPTION? 

Notation:  CD  =  coding  delay,  M  =  master,  S  =  slave,  LT  = 
latitude,  LN  =  longitude  and  BL  is  the  station  pair  baseline 
plus  the  secondary  phase  correction  in  microseconds. 
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4a.   Store  Data  in  Extended  Memory  -  XS 

With  the  extended  memory  module  in  the  HP-41CV,  station 
data  (in  R16  -  R25 ) ,  which  has  been  input  using  the  MI,  CR,  XR 
or  TR  routine,  can  be  transferred  to  the  module  using  the 
XEQ**XSn  command  or  by  pressing  the  appropriate  user  defined 
key.     Example  using  9940W. 

See  Key  in  Press 

1.  XEQnXS" 

2.  NEXT  OPTION? 

Should  the  station  pair  already  be  in  extended  memory,  the 
message  DUP  FL  (duplicate  file)  will  be  displayed.  If  needed, 
the  duplicate  file  may  be  erased  using  the  XD  routine. 

To  proceed  with  the  remaining  examples  it  is  recommended 
that  you  also  store  the  9940Y  station  pair. 

4b.   Recall  Data  from  Extended  Memory  -  XR 

With  the  extended  memory  module  installed,  XEQ"XRn  or 
press  the  appropriate  user  defined  key  to  input  the  data  for  a 
station  pair. 

Example  using  9940W.  (Note:  The  notation  <9>,  for  example, 
means  that  you  must  press  the  gold  shift  key  and  then  the  9 
key.) 

Se_e_  Key  in        £xe_s_a 

1.  XEQ"XRn 

2.  ID?  <9><9><4><0>W        R/S 

3.  NEXT  OPTION? 
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4c.   Delete  Data  from  Extended  Memory  -  XD 

XEQ"XD"  or  press  the  appropriate  user  defined  key  to 

delete  the  specific  station  pair  data  from  the  extended  memory 

module. 

Example  using  9940W. 

£££  Key  in        Zlzss. 

1.  XEQ"XDn 

2.  ID?  <9><9><4><0>W       R/S 

3.  NEXT  OPTION? 

The  message  FL  NOT  FOUND  will  be  displayed  if  the  file  you 
wish  to  delete  in  not  in  the  extended  memory. 

Note:  The  extended  functions/  memory  module  will 
accommodate  the  data  for  11  station  pairs.  The  extended  memory 
module  will  accommodate  the  data  for  an  additional  22  station 
pairs. 
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5a.   Store  Data  in  the  Tape  Cassette  -  TS 

With  the  tape  cassette  attached  to  the  HP-41CV,  station 
data  (in  R16  -  R25),  which  has  been  input  using  the  MI,  CR,  XR 
or  TR  routine,  can  be  transferred  to  the  tape  using  the  XEQnTS" 
command  or  by  pressing  the  appropriate  user  defined  key. 

Example  using  9940W. 

See  Key  in  Press 

1.  XEQnTS" 

2.  NEXT  OPTION? 

Should  the  station  pair  already  be  in  extended  memory,  the 
message  DUP  FL  NAME  (duplicate  file  name)  will  be  displayed.  If 
needed,  the  duplicate  file  may  be  erased  using  the  TD  routine. 

To  proceed  with  the  remaining  examples  it  is  recommended 
that  you  also  store  the  9940Y  station  pair. 

5b.   Recall  Data  from  the  Tape  Cassette  -TR 

With  the  tape  cassette  attached,  XEQ"TR"  or  press  the 

appropriate  user  defined  key  to  input  the  data  for  a  station 

pair . 

Example  using  9940W.  (Note:  The  notation  <9>,  for  example, 

means  that  you  must  press  the  gold  shift  key  and  then  the  9 

key.) 

See.  Key  in        elzss. 

1.  XEQnTR" 

2.  ID?  <9><9><4><0>W        R/S 

3.  NEXT  OPTION? 
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5c.   Delete  Data  from  the  Tape  Cassette  -  TD 

XEQ"TDn  or  press  the  appropriate  user  defined  key  to 

delete  the  specific  station  pair  data  from  the  tape  cassette 
Example  using  9940W. 

See.  Key  in        Pliers 

1.  XEQnTDn 

2.  ID?  <9><9><4><0>W       R/S 

3.  NEXT  OPTION? 

The  message  FL  NOT  FOUND  will  be  displayed  if  the  file  to  be 
deleted  is  not  on  the  tape. 
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6.    Loran-C  Fixing  Routines  FI  and  AS 

Given  the  indicated  time  delay  (ITD)  from  two  station 
pairs  which  form  a  triplet,  a  Loran-C  fix  is  obtained. 

Example:  Load  9940W  and  9940Y  into  the  calculator  using 
the  MI,  CR,  XR  or  TR  routine.  The  ITD  on  9940W  is  16019 
microseconds  and  the  ITD  of  9940Y  is  42585  microseconds.  Where 
are  you? 

see  Key  in        p_r_£££ 

1.  XEQ"FIn 

2.  ITD:  9940W  16019            R/S 

3.  ITD:  9940Y  42585            R/S 

4.  LAT:  39.1419  R/S 

5.  LON:  115.5052  R/S 

6.  NEXT  OPTION? 

Since  you  are  on  a  boat,  you  know  that  you  cannot  be  in 

central  Nevada  at  39dl4'19"  North  and  115d50,52"  West.  Every 

Loran-C  fix  has  two  solutions,  so  in  this  case  you  must  use  the 
alternate  solution. 

£e_£  Key  in        Rl&^s. 

7.  XEQ"AS" 

8.  LAT:  35.0001  R/S 

9.  LON:  125.0009  R/S 
10.  NEXT  OPTION? 

This  is  the  proper  solution  at  almost  exactly  35  degrees 
North  and  125  degrees  West.  Note  that  annunciator  3  (Flag  3) 
shows  in  the  display  indicating  the  alternate  solution.  If  you 
should  now  repeat  from  Step  1,  you  will  obtain  the  proper 
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solution  immediately. 

The  message  "E:  NO  TRIPLET"  will  appear  following  Step  1 
if  the  data  do  not  comprise  a  valid  triplet.  The  latitudes  and 
longitudes  of  each  station  pair  at  the  vertex  must  agree 
exactly.  Should  this  error  occur,  use  the  ED  routine  to  review 
the  resident  station  data. 

The  message  "E:  ITD  ERROR"  will  appear  following  Step  2  or 
3  indicating  that  the  ITD  you  keyed  in  is  inconsistent  with  the 
station  parameters.  Press  R/S  to  be  requeried  for  the  ITD. 
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7a.   Distance  and  Heading  Routines  DN  and  HP 

If  you  know  the  latitude  and  longitude  of  your 
destination,  you  may  key  these  in  and  then  see  how  far  your  fix 
is  from  your  destination  and  what  the  geodesic  heading 
(similiar  to  great  circle  heading)  is  to  your  destination. 

Example:  Your  destination  is  Moss  Landing  at  about  36d48'N 
and  121d47'W.  Your  current  fix  is  35dN  and  125dW  (see  the  FI-SA 
example) .  First,  key  in  your  destination. 

See  Key  in        P_r_es_s_ 

1.  XEQnDN" 

2.  DEST  LAT?  36.48  R/S 

3.  DEST  LON?  121.47  R/S 

4.  NEXT  OPTION?  R/S 

The  destination  is  now  stored  in  the  calculator  and  will  remain 
unchanged  until  you  use  either  the  DN  or  DH  options.  Also,  the 
latest  fix  is  stored  and  will  remain  unchanged  until  you  use 
either  the  FI,  AS  or  DH  options.  Now,  find  the  distance  and 
bearing  from  the  latest  fix  (see  the  FI-AS  example)  to  Moss 
Landing. 

s_££  Key  in  pjiess. 

1.  XEQnHD" 

2.  N.MI:    190.38  R/S 

3.  BRG:       54.3411  R/S 

4.  NEXT   OPTION?  R/S 

So  the  distance  to  Moss  Landing  is  190.38  nautical  miles 
at  a  heading  of  54d34,ll". 
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7b.   Distance  and  Heading  Routine  -  DH 

Given  the  latitude  and  longitude  of  an  origin  and 
destination,  this  routine  will  find  the  distance  and  heading 
from  one  to  the  other. 

Example:  How  far,  and  in  what  direction,  is  Corvallis, 
Oregon  (44d34'N,  123dl6 'W)  from  Cupertino,  California  (37dl9N, 
122d02'W)? 

£££  Key  in  Press 

1.  XEQ"DH" 

2.  ORIG  LAT?  37.19  R/S 

3.  ORIG  LON?  122.02  R/S 

4.  DEST  LAT?  44.34  R/S 

5.  DEST  LON?  123.16  R/S 

6.  N.MI:  438.32  R/S 

7.  BRG:   353.0259  R/S 

8.  NEXT  OPTION? 

Thus  the  distance  is  438.32  nautical  miles  and  the 
direction  of  Corvallis  from  Cupertino  is  353d02'59". 
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8.    ITD  Prediction  Routine  -  PR 

As  an  aid  to  identification,  this  routine  will  allow  the 
user  to  determine  what  ITD's  should  be  received  at  a  given 
location. 

Example:  Suppose  that  you  know  you  are  somewhere  near 
latitude  35  North  and  longitude  125  West  but  are  not  sure  what 
ITD's  you  should  be  receiving  from  9940W  and  9940Y.  To 
determine  these  ITD's,  proceed  as  follows: 
See 
1. 

2.  LAT? 

3.  LON? 

4.  9940W:  16019.35 

5.  9940Y:  42584.71 

6.  NEXT  OPTION? 

You  should  expect  to  receive  an  ITD  of  16019.35  from  9940W 
and  an  ITD  of  42584.71  from  9940Y. 


Key  in 

Press 

XEQnPR" 

35 

R/S 

125 

R/S 

R/S 

R/S 
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9.    Calibration  Routine  -  CA 

This  routine  will  allow  the  user  to  calibrate  the  Loran 
data  in  the  calculator  to  a  known  position  when  the  indicated 
time  delay  (ITD)  is  known  for  each  station  pair. 

Example:  Suppose  you  are  receiving  an  ITD  of  16308  from 
9940W  and  42800  from  9940Y.  These  ITD's  would  tell  you  that 
your  location  is  36d47'55"N  and  121d47 ' 11 "W.  However,  you  know 
that  your  position  is  bench  marked  to  be  at  36d47'36"N  and 
121d46 '58 "Wr  and  you  wish  to  calibrate  your  calculator  so  that 
the  ITD's  of  16308  and  42800  will  give  you  the  latter  fix 
instead  of  the  former.  Proceed  as  follows: 
See 
1. 

2.  LAT? 

3.  LON? 

4.  ITD  9940W: 

5.  ITD  9940Y: 

6.  NEXT  OPTION? 

Entering  16308  and  42800  into  the  FI  routine  will  now  give 
you  a  fix  at  36d47'35nN  and  121d46'59nW.  The  small  discrepancy 
between  this  fix  and  the  bench  mark  is  due  to  assumptions  made 
in  the  fixing  algorithm. 

Calibration  is  achieved  by  modifying  the  Master/Slave 
baseline  (BL  in  the  CE  and  ED  routines) .  See  Section  F. 


Key  in 

Press 

XEQnCA" 

36.4736 

R/S 

121.4658 

R/S 

16308 

R/S 

42800 

R/S 
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10.   Switch  Data  Registers  Routine  -  SW 

The  SW  utility  allows  the  user  to  swap  the  station  data 
stored  in  R16  -  R25  with  the  data  stored  in  R26  -  R35.  Whenever 
the  MI,  XR  or  TR  routine  is  used,  the  SW  routine  is  invoked 
prior  to  the  loading  of  the  data;  the  incoming  data  are  then 
placed  in  R16  -  R25 .  The  consequence  is  that  the  first  station 
pair  data  reside  in  R26  -  R35  and  the  second  station  pair  data 
reside  in  R16  -  R25 . 

One  user  application  of  SW  would  be  to  change  the  order  of 
the  station  ID  query  in  the  FI  routine  (this  also  affects  the 
order  of  determination  of  the  solution  and  alternate  solution) . 
Another  user  application  would  be  to  output  both  resident 
station  pairs  to  card,  extended  memory  or  tape  cassette  using 
the  CS,  XS  or  TS  routines,  respectively.  To  accomplish  this, 
first  use  the  CS,  XS,  or  TS  routine;  then  XEQnSW"  (note  that 
the  ID  of  the  station  data  in  R16  -  R25  appears  in  the  display 
instead  of  the  NEXT  OPTION?  prompt) ;  and  finally  use  the  CS,  XS 
or  TS  routine  once  more. 
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Ha.  Recording  the  Loran-C  Program  onto  Magnetic  Cards 

(1)  Attach  the  card  reader  to  the  HP-41CV. 

(2)  Place  the  calculator  in  the  USER  mode. 

(3)  Press  the  PRGM  key  to  place  the  calculator  in  the 
program  mode. 

(4)  Pass  one  side  of  a  blank  magnetic  card  through  the 
card  reader.  Then  follow  the  display  prompts  until 
nine  program  cards  (17  tracks)  have  been  recorded. 

(5)  Press  the  PRGM  key  once  more  to  leave  the  program 
mode. 

(6)  To  record  a  status  card,  XEQ"WSTS".  Then  pass  a  blank 
card  through  the  card  reader  following  the  display 
prompts. 

lib.  Recording  the  Loran-C  Program  onto  Magnetic  Tape 

(1)  Attach  the  tape  cassette  to  the  HP-41CV. 

(2)  Place  the  calculator  in  the  USER  mode. 

(3)  Press  the  alpha  key,  key  in  the  word  LORANC,  press  the 
alpha  key  once  more,  and  then  XEQnWRTP". 

(4)  To  record  the  program  status,  press  the  alpha  key,  key 
in  the  word  STATUS,  press  the  alpha  key  once  more,  and 
then  XEQ"WRTSn. 
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12a.  Loading  the  Loran-C  Program  from  Magnetic  Cards 

(1)  Attach  the  card  reader  to  the  HP-41CV. 

(2)  Clear  program  memory:  Turn  the  calculator  off,  then, 
while  pressing  the  left  arrow  (erase)  key  down,  turn 
the  calculator  on. 

(3)  Place  the  calculator  in  the  USER  mode. 

(4)  Read  in  the  STATUS  card.  The  status  card  will  set  the 
calculator  to  SIZE  42. 

(5)  Read  in  the  nine  program  cards  (17  tracks). 

12b.  Loading  the  Loran-C  program  from  Tape  Cassette 

(1)  Attach  the  tape  cassette  to  the  HP-41CV. 

(2)  Clear  program  memory.  (See  Step  2  above) . 

(3)  Place  the  calculator  in  the  USER  mode. 

(4)  Press  the  alpha  key,  key  in  the  word  STATUS,  press  the 
alpha  key  once  more,  and  then  XEQ"READSn.  The  status 
file  will  set  the  calculator  to  SIZE  42. 

(5)  Press  the  alpha  key,  key  in  the  word  LORANC,  press  the 
alpha  key  once  more,  and  then  XEQnREADP". 
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E.   The  Loran-C  Fixing  Algorithm 

The  principles  of  Loran  lines  of  position  (LOP's)  and 
fixing  are  adequately  covered  in  Reference  1  and  will  not  be 
repeated  here. 

The  basic  Loran-C  equation  [Ref.  4]  can  be  written  as 

ITD  =  [Ts  +  p(Ts)]  -  [TM  +  P(TM)]  +  [TB  +  p  (Tg)  ]  +  6         (1) 

where 

ITD  is  the  "indicated  time  difference"  in  microseconds, 
T   is  the  distance,  in  microseconds,  from  the  master 

to  the  receiver, 
Tq  is  the  distance,  in  microseconds,  from  the  slave  to 

the  receiver, 
TD  is  the  distance,  in  microseconds,  between  the  master 
and  the  slave, 
6  is  the  assigned  station  pair  coding  delay,  in  micro- 
seconds, and 
p(T)  is  the  secondary  phase  correction,  in  microseconds, 
for  a  surface  seawater  path  of  length  T  . 
The  quantity 

At  =   TB  +  p(TB)   +  6 

is  a  constant  for  each  master/slave  pair.  The  quantity   T    is 

computed  from  the  positions  of  the  master  and  slave  using  the 

reverse  solution  algorithm  (Section  H)  at  the  time  of  manual 
data  input  (Routine  MI) . 
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The  following  World  Geodetic  System  1972  (WGS  72)  values 
have  been  adopted  for  Loran-C  navigation  [Ref.  4]: 

vn  =  299792458  meters/second  is  the  velocity  of  light 

in  free  space, 

n  =  1.000338  is  the  index  of  refraction  of  the 

surface  of  the  earth  for  standard  atmosphere 

and  100  kHz  electromagnetic  waves, 

a   =  6378135.000  meters  is  the  equatorial  radius  of 
e  ^ 

the  earth 
and        f  =  1/298.26  is  the  flattening  factor  (1-b/a  ,  where 

b   is  the  polar  radius)  of  the  earth. 

Accurate  formulas  for  computing  the  secondary  phase  cor- 
rection p(T)  are  contained  in  Reference  4,  but  for  use  in  the 
HP- 4 1CV,  the  form 


p(T)  =  aQ/T  +  ax  +  a2T  (2) 


is  used,  where   T   is  in  microseconds  and 

1.  For   T  >  537  ysec: 

aQ  =  129, 

a,  =  -0.408,  and 

a„  =  0.0006458. 

2.  For   T  <  537  psec: 

aQ  =  2.74, 

a   =  -0.011,  and 

a2  =  0.00033. 
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On  the  surface  of  a  sphere,  a  hyperbolic  line  of  position 
(LOP)  can  be  represented  by  the  equation  [Ref.  1,  page  175] 

4-^„  v.  _    cos  2a  -  cos  2c  ,_. 

tan  r  =  — : ~ : ^ —  (3) 

sin  2c  cos  co  +  c,    sin  2a 

where  the  origin  of  the  coordinate  system  is  at  the  prime  focus 
of  the  spherical  hyperbola,  2c  is  the  spherical  arc  joining  the 
foci,  2a  is  a  constant  for  any  one  hyperbola,  and   r   and   co 
are  the  spherical  coordinates  of  a  point  on  the  hyperbola.   If 
the  base  line  of  the  coordinate  system  is  the  arc  joining  the 
foci  the   oo   is  the  spherical  polar  angle  from  the  baseline  to 
a  point   P   on  the  spherical  hyperbola  and   r   is  spherical 
polar  distance  (or  arc)  from  the  prime  focus  to   P  .   Using  the 
Loran  system  we  take   C  =  +1   if  the  prime  focus  is  at  a  master 
station  and   £  =  -1   if  the  prime  focus  is  at  a  slave  station. 

If  we  let   v  =  v^/n   be  the  velocity  of  100kHz  electro- 
magnetic radiation  at  the  earth's  surface  then,  for  a  spherical 
earth,  we  can  relate  the  parameters  in  Equations  1  and  3  as 
follows : 


and 


2c  =  vTB/ae  , 


2a  -  v(T  -  TM)/ae 


Using  the  spherical  approximation  for  now,  we  see  that  2c  is 
known  for  any  Loran  pair.   The  "indicated  time  delay"  ITD  is 
measured  by  the  receiver  at  point   P  ,  and  to  determine  a 
hyperbolic  line  of  position  we  must  determine  2a,  but  T   -  T 
cannot  be  computed  from  Equations  1  and  2.   If  a„  were  zero 
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in  Equation  2,  then  it  would  be  possible  to  determine  T   -  TM 
uniquely.   As  an  approximation  we  use  the  following  parameters 
in  Equation  2 : 

aQ  =  0  , 
a,  =  -0.321  , 
and  a2  =  0.000635  . 

These  values  have  been  obtained  by  setting  a„  =  0  and  determin- 
ing a,  and  a~  by  linear  regression  of  the  T  >  537  values  over 
the  interval  of  1000  <  T  <  8000.   This  approximation  is  quite 
good  (within  0.03  us)  for  distances  up  to  10,000  microseconds 
where  small  changes  in  the  LOP ' s  can  cause  large  position  errors. 
At  short  distances  the  error  increases  from  0.05  us  at  1000  us 
to  0.58  ys  at  10  us;  although  these  errors  are  large  for  small 
distances,  the  LOP ' s  are  not  as  sensitive  to  these  changes  as 
they  would  be  at  large  distances.   When  this  approximation  is 
substituted  into  Equation  1,  we  obtain 

[Ts  +  ax  +  a2Tg]  -  [TM  +  a1  +  a^]  =  ITD  -  At  , 
or 

Ts  -  TM  =  (ITD  -  At)/(1  +  a2)  (4) 

and  hence  2a  =  v(T  -T  )/a   is  determined  for  use  in  the 
spherical  approximation. 

Consider  a  Loran-C  triplet  with  the  master  stations  co- 
located.   Let   £,   and   ^„   denote  the  azimuth  angles  of  slave 
1  (S,  )  and  slave  2  (S2) ,  respectively,  measured  from  North 
toward  the  East  from  the  master  stations  (M)  (see  Figure  3) . 
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Further,  let   a   and   r   be  the  azimuth  and  spherical  polar 

arc  (distance)  of  the  receiver  (R)  from  M  .   For  this  geometry, 

Equation  3  can  be  written  in  the  form 


B. 

tan  r.  =  - -. — = — r — ■ — ,  (5) 

1    C .  cos (a  -  £ . )  +  A.   '  K    ' 

i  i     l 

where 

A.  =  C •  sin  2a .  , 
11        l  ' 

B.  =  cos  2a.  -  cos  2c  , 
l  l  l  ' 


and 


C.  =  sin  2c. 
l  l 


•  th 


for  the   i     Loran  pair,  i  =  1,2.   Since  r,  =  r„  =  r,  we  can 
eliminate   tan  r   between  the  two  equations.   The  resulting 
equation  can  be  rewritten  as 


where 


and 


C  cos  a  +  S  sin  a  =  K  ,  (6) 


C  =  B,C2  cos  ^2  -  B2C,  cos  £,  , 


S  =  B,C2  sin  E,2    -   B2C-,  sin  £,  , 


K  =  B2A1  -  BXA2  . 


If  we  define   p  >  0   and  y      by  the  equations 

p  cos  y  =  C  ,  (7) 

and  p  sin  y  =   S  , 

then 


=  Jc2 


p  =  /c  +  s 

and  Y  =  qatn  (S,C) 
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Figure  3.   Geometry  of  a  Loran  Triplet  and  a  Receiver 
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Here  the  function  qatn(y,x)  is  the  arctangent  of  y/x  adjusted 
for  the  proper  quadrant  according  to  the  signs  of   x   and   y  . 
A  compact  form  of  this  function  is 


where 
and 


qatn(y,x)  =  tan 

t(z)  =  1 
t(z)  =  0 


-1 


x  +  10  9  t(x  =  0?) 


+  TTt  (X  <  0?) 


when   z   is  true 
when   z   is  false. 


When  convenient  we  will  use  the  notation  qatn(y/x)  inter- 
changeably with  qatn(y,x).   Now  we  can  substitute  Eq.  (7)  into 
Eq.  (6)  and  solve  for 

a  =  y  ±  cos   (K/p) 


(8) 


to  obtain  the  azimuth  angle   a   of  the  two  points  of  inter- 
section of  the  spherical  hyperbolic  LOP's.   Finally  we  can 
obtain  a  value  for   r   by  substituting  each   a   into  Eq.  (5) 
for  either  i  =  1  or  i  =  2 .   We  find  that 


r  =  qatn 


B 


C.  cos (a-£ . )  +  A. 
l         l      l 


for  i  =  1  or  2 


(9) 


The  distance  and  azimuth  from  M  or  the  triplet  vertex  can  be 
converted  into  the  latitude  and  longitude  of  the  two  possible 
positions  of   r  . 

The  fixing  algorithm  then  uses   a   and   r   in  the  direct 
solution  algorithm  of  spheroidal  geodesy  (Section  G) . 
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F.   The  Calibration  and  ITD  Prediction  Algorithms 

Calibration  can  be  achieved  when  an  ITD  is  measured  at  a 
known  bench  marked  position.   From  the  bench  marked  position 
and  the  known  master  and  slave  positions,  the  quantities 

T„  +  p(T„)   and   T   +  p(T  )   can  be  computed  using  the  reverse 

M   ^   M         s    r   s  r 

solution  algorithm  (Section  H)  and  the  accurate  secondary  phase 

correction  formula  (Eq.  2) .   Equation  1  can  then  be  solved  for 

T_,  +  p(TD)   to  obtain  a  modified  baseline.   This  modified  base- 
B        B 

line  is  stored  and  then  used  instead  of  the  true  baseline  in 
subsequent  computations.   The  affect  on  the  accuracy  of  fixes 
using  this  modified  baseline  with  positions  far  removed  from 
the  bench  mark  has  not  studied. 

The  ITD  prediction  algorithm  is  a  direct  application  of 
Equation  1.   A  known  position,  together  with  the  known  master 
and  slave  positions,  is  used  to  compute  the  quantities 
T   +  P(TM)   and   T   +  p(T  ).   When  these  values,  along  with  a 
computed  or  calibrated  baseline,  TD  +  p(TD)  ,  are  substituted 

B         B 

into  Equation  1,  a  predicted  ITD  is  obtained. 
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G.   The  Direct  Solution  Algorithm 

This  direct  solution  algorithm  is  a  modification  of  the 
second  order  in  flattening  (f)  algorithm  given  by  Thomas 
[Ref.  5,  pp.  7-8].   Thomas'  notation  has  been  followed  as 
closely  as  possible  for  ease  of  comparison  of  the  algorithms. 
The  qatn  function  is  defined  in  a  previous  section.   East 
longitudes  and  South  latitudes  are  negative.   We  are  given  the 
point  P,  (c|),,A,)  on  the  spheroid,  where  $,  ,  A,  are  the  geodetic 
latitude  and  longitude  (geographic  coordinates) ;  the  forward 
azimuth  a,  ~  and  distance   S   to  a  second  point  P_  (<j>~  ,  A~)  ;  and 
from  these  we  are  to  find  the  geographic  coordinates  <£,,,  A~  and 
the  back  azimuth  a~,  .   The  given  quantities  are   <}>,  ,  A,,  a,  ~ 
and   S  .   No  assumptions  about  the  relative  location  of   P, 
and   P~   are  required.   The  modified  direct  solution  algorithm 
is : 

8.  =  tan   [(l-f)tan  6,],  M  =  -  sin  a,-  cos  6,  , 
C1  =  fM,  C2  =  f(l-M2)/4  , 

D  =  d-C2) (1-C2  -  C±M) ,  P  =  C2El  +  (1/2)C1M]/D  , 

N  =  cos  6,  cos  a. „  ,  a,  =  qatn(N,  sin  6,)  , 

d  =  S/(a  D)  ,  u  =  2(a1-d)  ,  W  =  1-2P  cos  u  , 

2  2 

V  =  cos(u+d)  ,  X  =  C2  sin  d  cos  d(2V  -1)  , 

Y  =  2PVW  sin  d,Aa=d+X-Y, 

2  2  1/2 

K  =  [M   +  (N  cos-Aa  sin  6,  sin  Aa)  ]     , 

6-  =  tan"  [(sin  6,  cos  Ao  +  N  sin  Aa)/K]  , 

An  =  qatn(-  sin  Aa  sin  a,2  ,  cos  6,  cos  Aa  -  sin  6,  sin  Aa  cos  a-i2)  , 

H  =  C,  (1-C2)  Aa  -  C1C2  sin  Aa  cos^a^^  -  Aa)  , 

A2  =  A,  +  An  -  H  , 
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'    =  qatn[-M,  -  (N  cos  Aa  -  sin  0.  sin  Aa) ]  , 


<p2   =    tan  ""-[tan  62/(l-f)] 

Details  of  the  modifications  made  to  Thomas'  algorithm  are 
contained  in  Reference  3.   The  algorithm  above  has  been  further 
modified  so  that  Eastern  longitudes,  rather  than  Western 
longitudes,  are  negative. 
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H.   The  Reverse  Solution  Algorithm 

This  reverse  solution  algorithm  is  a  modification  of  the 
second  order  in  flattening  (f)  algorithm  given  by  Thomas 
[Ref.  5,  pp.  8-10].   Thomas'  notation  has  been  followed  as 
closely  as  possible  for  ease  of  comparison  of  the  algorithms. 
The  qatn  function  is  defined  in  a  previous  section.   East 
longitudes  (A)  and  South  latitudes  (<J>)  are  negative.   We  are 
given  the  points   P,  (<}>,,  A.),  P2(t}>2,A2)  on  the  spheroid  and  are 
to  find  the  distance   S   between  the  points  and  the  forward 
and  back  azimuths,  a,„  and  a21  •   Given  quantities  are 
<J>i  /  ^-i  i    $?    an<3  ^->  •   No  assumptions  about  the  relative  location 
of  P-,  and  P~  are  required.   The  modified  reverse  solution 
algorithm  is: 

6.  =  tan_1[  (l-f)tan  (Jk]  ,  i  =  1,2  , 

AA  =  A2  -  A1  ,  AG^  =  (92  -  6^/2  ,  6m  =  (Q±   +    QJ/2    , 

H  =  cos2  A6   -  sin2  6   ,  L  =  sin2  A6   +  H  sin2 (AX/2)  , 

mm  m 

_  1  2        2 

d  =  2  sin   /L  ,  U  =  2  sin   6   cos   A0  /(1-L)  , 

m        m 

V  =  2  sin2  A6   cos2  0/L,X=U+V,Y=U-V, 
mm 

2 
T  =  d/sin  d  ,  D  =  4T   ,  E  =  2  cos  d  ,  A  =  DE  , 

C  =  T  -  (A-E)/2,  n   =  X(A+CX)  , 

B  =  2D,  n2  -  Y(B+EY)  ,  n3  =  DXY  , 

62d  =  f2(n1  -  n2  +  n3)/64  ,  6±d   =  f(TX-Y)/4  , 

S/a   =  (T  -  6?d  +  60d)  sin  d  ,  M  =  32T  -  (20T-A)X  -  (B+4)Y  , 
e         1      ^ 

F  =  2Y  -  E(4-X)  ,  G  =  fT/2  +  f2M/64  ,  Q  =  -  (FG  tan  AA)/4  , 

AA  '   =  (AA  +  Q)/2  , 
m 
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tn  =  qatn(sin  A0   cos  AA'  ,  cos  0   sin  AA')  , 

1  ^  ra       m        m       m 

t~  =  qatn(-  cos  A6   cos  AA1  ,  sin  9   sin  AA')  , 

2  ^  m       m        m       m 

a12  =  tl   +   t2    '    a21  =  ""^l  ~  fc2 


Details  of  the  modifications  made  to  Thomas'  algorithm 
are  contained  in  Reference  3.   The  algorithm  above  has  been 
further  modified  so  that  Eastern  longitudes,  rather  than 
Western  longitudes,  are  negative. 
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I.  Program  Accuracy 

The  direct  and  reverse  solution  algorithms  are  equivalent  to 
the  second  order  flattening  algorithms  given  by  Thomas  (Ref.  5); 
the  parameters  of  the  WGS  1972  spheroid  are  used.  The  reverse 
solution  algorithm  reproduces  the  baselines  provided  by  the 
Defense  Mapping  Agency  for  all  40  Loran-C  stations  to  within  0.15 
meters  (the  average  deviation  is  -0.031  meters,  DMA  minus  HP-41CV, 
with  a  standard  deviation  of  0.037  meters)  and  to  within  0.01 
microseconds,  including  the  secondary  phase  correction  for  an  all 
seawater  path.  The  reverse  solution  algorithm  is  also  used  to 
generate  predicted  ITD's;  these  are  presumed  to  be  within  0.01 
microseconds,  also. 

The  fixing  algorithm  uses  the  direct  solution  algorithm  with 
the  azimuth  and  distance  of  the  fix  from  the  vertex  of  the  Loran-C 
triplet  computed  from  Equations  8  and  9  as  inputs.  Equations  8  and 
9  are  based  upon  spherical  geometry  and  include  an  approximation 
to  the  secondary  phase  correction  for  an  all  seawater  path.  The 
largest  source  of  error  is  the  assumption  that  the  azimuth  and 
distance  to  the  fix  are  accurately  represented  by  this  spherical 
approximation.  This  approximation  has  not  been  rigorously  tested, 
however  it  is  possible  to  use  the  reverse  solution  algorithm  to 
predict  the  ITD's  that  will  be  received  at  a  given  position  and 
then  to  enter  these  ITD's  into  the  fixing  algorithm  to  determine 
how  accurately  the  fixing  algorithm  reproduces  the  original 
position.  The  distance  between  the  fix  and  the  original  position 
can  be  determined  using  the  HD  algorithm.  Tables  1,  2  and  3  were 
producted  in  this  manner.  Similiar  tables  with  different  station 
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pairs  are  given  in  References  3,  6  and  7.  It  is  felt  that  the 
results  are  accurate  enough  to  not  warrant  the  inclusion  of  an 
iterative  improvement  routine.  Of  the  samples  in  Tables  1,  2  and 
3,  the  largest  error,  0.22  n.mi.,  is  the  first  entry  in  Table  1. 
From  the  chart  LCNC-2,  it  estimated  that  the  angle  of  intersection 
of  the  two  hyperbolic  lines  of  position  is  about  5  degrees  and  so 
an  error  of  0.22  n.mi.  should  not  be  unexpected. 
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Table  1.   Colocated  Master  Stations 


Position 

Predicted  ITD's 

HP- 4 1CV  Fix 

Error 

Lat 

Long 

9940W 

9940X 

Lat(N) 

Long(W) 

n.mi 

31° 

123° 

16413.28 

27570.93 

30o59'47" 

123°00'03" 

0.22 

37° 

126° 

15610.11 

27020.50 

36°59'53" 

126°00'09" 

0.17 

42° 

129° 

13881.78 

27285.58 

42o00f00" 

129o00'00M 

0 

44° 

132° 

13180.89 

27371.19 

44°00,00" 

132°00'01" 

0.01 

48° 

135° 

12301.25 

27552.06 

48°00*01" 

135°00'03" 

0.03 

50° 

138° 

12068.67 

27584.22 

50°00f01M 

138o00'04M 

0.05 

Table  2.   Colocated  Slave  Stations 


Position 

Predicted  ITD's 

HP-41CV  Fix 

Error 

Lat 

Long 

9940W 

5990Y 

Lat(N) 

Long(W) 

n.mi 

31° 

123° 

16413.28 

27177.18 

31°00'03" 

122°59'58n 

0.06 

37° 

126° 

15610.11 

27403.20 

37°00»01" 

125°59'59" 

0.01 

42° 

129° 

13881.78 

27955.45 

42°00,00M 

129o00'00M 

0 

44° 

132° 

13180.89 

28512.90 

44°00'00" 

132o00'00M 

0 

48° 

135° 

12301.25 

29413.61 

48°00,00" 

134°59'59" 

0.01 

50° 

138° 

12068.67 

29816.84 

50°00,00" 

137°59'59" 

0.02 
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Table  3.   Colocated  Master  and  Slave  Stations 


Position 

Predicted  ITD's 

HP-41CV  Fix 

Error 

Lat 

Long 

5930Y 

9960W 

Lat(N) 

Long(W) 

n.mi 

44° 

63° 

29864.46 

11685.15 

44°00,00" 

63°00,00" 

0 

41° 

66° 

30585.61 

12946.91 

41o00'00" 

66°00'00" 

0 

39° 

69° 

31020.46 

14111.31 

39°00f00" 

69o00'00" 

0 

35° 

72° 

31064.57 

15139.48 

35°00'00M 

72°00'00M 

0 

30° 

75° 

31040.82 

15610.46 

29°59'59" 

75°00,01" 

0.02 

26° 

78° 

31106.20 

15858.46 

25°59'57" 

78o00,02'* 

0.05 
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Appendix  A:   Program  Storage  Allocations ,  Flag  Usage 
and  Program  Listing 

Program  Size:   4  2 

Registers: 

R00  -  R13:   Scratch  storage 

R14:   Flattening 


2nd  Sta. 

Pair 

Variable 

R15: 

2a 

R16: 

ID 

R17: 

TB  +  p(TB) 

R18: 

2c 

R19: 

CD 

R20: 

0  master 

R21: 

A  master 

R22: 

a  master-slave 

R23: 

G  slave 

R24: 

A  slave 

R25: 

a  slave-master 
2a 

1st  Sta.  Pair 


R26 
R27 
R28 
R29 
R30 
R31 
R32 
R33 
R34 
R35 
R36 


R37 
R38 
R39 
R4  0 
R41 


6  fix 

A  fix 

a  fix-dest 

0  dest 

A  dest 
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Flag  Usage: 

FOO:   1.  Requery  erroneous  ITD  at  LBL06 

2.  Echo  two  data  sets  at  LBL  ED 

F01:  PR  and  CA  interlock 

F03:  Set  for  2nd  solution 

F05:  Set  if  slave  at  vertex  of  2nd  station  pair 

F06:  Set  if  slave  at  vertex  of  1st  station  pair 

F0  7:  Set  for  XMEM  functions,  clear  for  tape  functions 

F08:  DN  and  DH  interlock 

F09:  Set  if  DMS  conversion  required 

F10:   1.  SW  and  AS  interlock 

2.  Input  interlock  in  FI   vertex  check 

F14:  PR  and  CA  loop  control 
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Appendix  B;     station  coverage 


Station    No.  of  Pairs 

Location 

4990 

2 

Central  Pacific 

5930 

2 

East  Coast,  Canada 

5990 

3 

West  Coast,  Canada 

7930 

3 

North  Atlantic 

7960 

2 

Gulf  of  Alaska 

7970 

4 

Norwegian  Sea 

7980 

4 

Southeast  U.S.A. 

7990 

3 

Mediterranean  Sea 

8970 

3 

Great  Lakes 

9940 

3 

West  coast,  U.S.A. 

9960 

4 

Northeast  U.S.A. 

9970 

4 

Northwest  Pacific 

9990 

3 

North  Pacific 
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Appendix  C:   LORAN-C  STATION  DATA 

The  following  list  contains  the  pertinent  parameters  for  each 
Loran-C  station  pair.  This  list  was  compiled  from  the  data  in 
Reference  4.  Each  column  contains  the  following  information: 

1.  The  Loran-C  station  pair  designator. 

2.  The  coding  delay. 

3.  The  master  station  latitude. 

4.  The  master  station  longitude. 

5.  The  slave  station  latitude. 

6.  The  slave  station  longitude. 

In  this  list,  positive  longitudes  are  West,  negative 
longitudes  are  East,  positive  latitudes  are  North  and  negative 
latitudes  are  South.  In  columns  3  through  6  the  latitudes  and 
longitudes  appear  to  be  in  decimal  form,  but  the  actual  format  is 
DDD.MMSSFF  (which  is  compatible  with  the  HP-41CV  D.MS  or  H.MS 
input  mode)  where 

DDD  designates  degrees, 
MM  designates  minutes, 
SS   designates  seconds,  and 
FF  designates  hundredths  of  seconds. 
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ID  CD  M£_  SS  LAT  SS  LQN 


4990X 

11000 

155.530970 

4990Y 

29000 

X  0 

- 

178.173020 

5930X 

11000 

46 

069.583909 

5930Y 

25000 

46 

053.102816 

5990X 

11000 

51 

131.151965 

5990Y 

27000 

51 

119.443953 

5990Z 

41000 

51 

127.212935 

7930W 

11000 

023.552175 

7930X 

Old  r"   ' 

007.042671 

7930Z 

053.102816 

7960X 

152.221122 

7960Y 

131.151965 

7970W 

O  C  f 

-008.173633 

7970X 

-014.274700 

7970Y 

460 

.    . 

+023.552175 

7970Z 

60  ' 

261 

+008.435869 

7980W 

11000 

302 

090.494360 

7980X 

23000 

**0I 

~  C  «  "I 

097.500009 

7980Y 

43000 

■  .  ■" 

-' 

080.065352 

7980Z 

59000 

30 

077.544676 

7990X 

11000 

38. 

-012.312996 

7990Y 

29000 

38. 

-027.520152 

7990Z 

47000 

-003.121590 

8970W 

11000 

39. 

085.100930 

8970X 

28000 

39, 

j  [Jq[J 

076.493386 

8970Y 

44000 

094.331847 

9940W 

11000 

39. 

799 

119.443953 

9940X 

27000 

.465699 

122.294453 

9940Y 

40000 

114.481743 

9960W 

11000 

A  *) 

067.553771 

9960X 

25002 

069.583909 

9960Y 

"3  Q 

077.544676 

9960Z 

087.291214 

9970W 

-153.58515 

9970X 

-143.430906 

9970Y 

>499 

-128.085621 

9970Z 

24 

-138.095523 

9990X 

11000 

3.105231 

9990Y 

29000 

+166 . 531447 

9990Z 

43000 

57 

,05 

57 

.262021 

+152.221122 

Appendix  D;   CQLQCATED  STATIONS 


5930  MASTER  -  9960W 
5930X  -  9960X 
5930Y  -  7930Z 
5990X  -  7960Y 
5990Y  -  9940W 
7930W  -  7970Y 
7930X  -  7970  MASTER 
7960X  -  9990Z 
7980Z  -  9960Y 
7980  MASTER  -  8970W 
8970X  -  9960  MASTER 
8970  MASTER  -  9960Z 


Slave  stations  are  denoted  with  a  letter 
suffix.  Master  stations  are  so  designated 
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Naval  Postgraduate  School 
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Office  of  Research  Administration  1 
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Naval  Postgraduate  School 

Monterey,  CA   93940 
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Naval  Postgraduate  School 
Monterey,  CA   93940 

Office  of  Naval  Research  2 

Fleet  Activity  Support  Division 
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800  North  Quincy  Street 
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Navy  Tactical  Support  Activiity  2 

P.O.  Box  1042 

Silver  Springs,  MD   20910 

Attn:  Mary  Ellen  Lannon 

COMPATWINGSPAC  2 

Naval  Air  Station 
Moffett  Field,  CA   94035 
Attn:  Code  51  and  Code  532 

COMPATWINGSLANT  2 

Naval  Air  Station 
Brunswick,  ME   04011 
Attn:  Code  N7 

Prof.  R.H.  Shudde,  Code  55Su  100 

Naval  Postgraduate  School 
Monterey,  CA  93940 
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